Momentum-space anisotropy and pseudogaps: a comparative cluster dynamical 
mean-field analysis of the doping-driven metal-insulator transition in the 

two-dimensional Hubbard model 
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Cluster dynamical mean-field calculations based on 2-, 4-, 8- and 16-site clusters are used to 
analyze the doping-driven metal-insulator transition in the two-dimensional Hubbard model. Com- 
parison of results obtained on different clusters enables a determination of those aspects of the 
physics that are common to all clusters and permits identification of artifacts associated with par- 
ticular cluster geometries. A modest particle-hole asymmetry in the underlying band structure is 
shown to lead to qualitatively different behavior on the hole-doped side than on the electron-doped 
side. For particle-hole asymmetry of the sign and magnitude appropriate to high-T c cuprates, 
the approach to the insulator from the hole-doping side is found to proceed in two stages from a 
high-doping region where the properties are those of a Fermi liquid with moderately renormalized 
parameters and very weak momentum dependence. As doping is reduced the system first enters 
an intermediate doping regime where the Fermi-liquid renormalizations are larger and the electron 
self-energy varies significantly around the Fermi surface and then passes to a small doping regime 
characterized by a gap on some parts of the Fermi surface but gapless behavior in other parts. On 
the electron-doped side the partially gapped regime does not occur, and the momentum dependence 
of the electron self-energy is less pronounced. Implications for the high-T c cuprates and for the use 
of cluster dynamical mean-field methods in wider classes of problems are discussed. 

PACS numbers: 71.10.Fd, 71.27.+a,71.30.+h,74.72.-h 



I. INTRODUCTION 

The evolution of a Mott (correlation-driven) insulating 
state as charge carriers are induced by doping is one of the 
basic questions in condensed-matter physics. It is rele- 
vant to the behavior of wide classes of materials^ and is of 
particular importance in the context of high-temperature 
copper-oxide superconductivity, where following an early 
insight of Anderson^ the materials are generally accepted 
to be doped Mott insulators, with the combination of 
"Mott" correlations and (quasi) two-dimensionality giv- 
ing rise to the novel properties. 

A crucial insight into the physics of doped Mott insu- 
lators was provided by Brinkman and Rice^ who argued 
that the low-temperature properties of the doped Mott 
insulator were in effect those of a Fermi liquid with a 
quasiparticle mass which diverged as the Mott insulating 
state was approached. The Brinkman-Rice picture pro- 
vides a reasonable description of data in many materials 1 - 
and has been refined theoretically over many years. In 
particular the development of the single-site dynamical 
mean-field theory 4 - provided a precise theoretical frame- 
work within which these results can be derived. The key 
physical assumptions required for Brinkman-Rice behav- 
ior are now understood to be the presence of strong corre- 
lations and the locality of correlation effects. The mathe- 
matical expression of the latter assumption is the momen- 
tum independence of the electron self-energy. The local- 
ity assumption becomes strictly valid for classes of lattice 



models in a limit of infinite lattice coordination number— 
and appears to provide a reasonable approximation to the 
behavior of many three-dimensional materials^ However, 
the momentum-independent self-energy approximation is 
likely to be less accurate for materials, such as high-T c 
superconductors, where the electronic properties are two- 
dimensional. 

Indeed, it was recognized soon after the discovery 
of high-T c superconductivity in the Cu02 perovskites 
that the materials required a description which went 
beyond the Brinkman-Rice/single-site dynamical mean- 
field theory approach. The evidence has become stronger 
over the years and it is beyond the scope of this pa- 
per to review the large experimental literature demon- 
strating this point. We do recall here the results of 
angle-resolved photoemission and related transport ex- 
periments which clearly illustrate this point. Studies 
of hole-doped cuprate materials with very high dopings 
(above the doping which maximizes the superconducting 
transition temperature) indicate a quasiparticle lifetime 
and velocity renormalization which are nearly isotropic 
around the Fermi surface^ consistent with inferences 
from magnetotransportii^ As the doping is reduced to- 
ward the insulating phase the behavior changes. Pho- 
toemission measurements on materials with "optimal" 
doping levels (near the carrier concentration which maxi- 
mizes the superconducting transition temperature) reveal 
an electron lifetime which varies dramatically around the 
Fermi surface^ again consistent with inferences drawn 



from interpretations of the interplane conductivity 1 ^ and 
from in-plane magnetoresistance measurements^ Mea- 
surements on "undcrdoped" cuprates (those with carrier 
concentrations even closer to the insulating phase) re- 
veal a "pseudogap" (reduction in electronic density of 
states) for momenta near the (0,ir)/(ir,0) points of the 
Brillouin zone but no pseudogap for states near the zone- 
diagonal (near (ir/2,ir/2))J^-~— Thus optimally doped 
and undcrdoped cuprate materials exhibit a self-energy 
with a strong momentum dependence, inconsistent with 
the Brinkman-Rice/single-site dynamical mean-field ap- 
proach. 

On the theoretical side the importance of going beyond 
the Brinkman-Rice description was also recognized early 
on. For example, the various forms of resonating valence 
bond ( "RVB" ) theories take intersite correlations explic- 
itly into account by expansion around a specific mean- 
field approximation and were shown to lead to a non- 
Brinkman-Rice doping dependence of the electron effec- 
tive mass and other Fermi liquid parameter o 14 ' 15 and to 
a pseudogap^ (for a recent review see Ref. Il7l ). Other 
authors have addressed the issue in the context of ana- 
lytical calculations based on the assumed importance of 
antiferromagnetioi^— or charge-density wave 2 ^ correla- 
tions. These, and other semi-analytic calculations, while 
demonstrating the importance of antifcrromagnctic cor- 
relations for the electron-doped cuprates, have not led 
to a consensus regarding the physics of the hole-doped 
cuprates, in part because they are based on approxima- 
tions which are uncontrolled (or are controlled in limits 
which are not clearly relevant to the actual materials) 
and more importantly because they are based on assump- 
tions about which correlations are physically relevant and 
which may be neglected. The appropriateness of the un- 
derlying assumptions about which physics to include have 
been the subject of debate. 

Over the last decade the development of "cluster" dy- 
namical mean-field methods^ 3 - - — has opened up a very 
promising new line of attack on the problem. These 
methods obtain an approximate solution of the full many- 
body problem in terms of the solution of an auxiliary 
N-site quantum impurity model coupled with a self- 
consistency condition. For a review, see Ref. [27| and 
Refs. 28.29. N = 1 corresponds to single-site dynamical 
mean-field theory with a momentum-independent self- 
energy; clusters of size N > 1 allow for some momentum 
dependence of the self-energy and thus enable the study 
of deviations from Brinkman-Rice behavior. As N — > oo 
one recovers the full model; however, the computational 
expense rises rapidly as the interaction strength and clus- 
ter size increase. An advantage of the methods is that 
no explicit assumption is made about the important of 
one kind of electronic correlation (spin density, charge 
density RVB) relative to another, but the possibility of 
potential biases associated with choice of cluster is an 
important issue which this paper aims to address. 

Important results obtained by cluster dynamical 
mean-field methods have included the demonstration 
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FIG. 1: Qualitative sketch of doping regimes for parameters 
considered in this paper. 



that in an appropriate doping and interaction range 
the two-dimensional Hubbard model can exhibit a 
pseudogap ) 30 ' 31 "Fermi arcs"^2r— and a variation in 
electronic properties around the Fermi surface ("nodal- 
antinodal differentiation" ) J^-33^1-39 rp ne me thods have 
been shown to yield a multistage approach to the Mott 
transition with the insulating phase being separated 
from the weakly correlated metallic phase by a "sector- 
selective" phase where some regions of the Brillouin zone 
are gapped and others are no t 37 ' 38 ' 40 ' 41 and a number of 
physical properties were shown to be in good agreement 
with experiment . 38 ' 39 ' 42 

Despite these successes, uncertainties remain. Most 
cluster dynamical mean-field papers analyze one specific 
cluster. Comparison of results obtained on clusters of 
different sizes has been undertaken only in a few special 



cases 



43-47 



mostly not directly relevant to the question 
of the doping-driven Mott transition. It has therefore 
not been clear which results are due to specific proper- 
ties of clusters and which results are representative of the 
physics of the full model. More generally, cluster dynam- 
ical mean-field calculations test the limits of present day 
computational abilities so that compromises are required 
between cluster size and the ranges of temperature, inter- 
action strength and carrier concentration to be studied. 
Little information is available in the literature to guide 
the choices which must be made. 

This paper has two main goals: to clarify the physics of 
the doping-driven Mott transition in two dimensions by 
identifying the robust physical features which can now be 
considered as established from cluster dynamical mean- 
field theory and to identify the differences between dif- 
ferent cluster sizes and geometries. We present a global 
examination of the one-electron properties of the doped 
Mott insulator, using clusters of all feasible sizes from 
2 up to 16 sites, respecting the lattice symmetry. The 
investigation is made possible by progress in algorithms, 
which have made the computations much more efficient, 
enabling surveys of wide ranges of parameter space for 
many clusters . 48 ' 49 We describe the electronic properties 
in detail and show that the phase diagram and physi- 
cal properties are to a surprising extent robust against 
choice of cluster size and geometry. Where differences 
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FIG. 2: (Color online) Momentum-space tiling used to define cluster approximations studied here: 2-site (leftmost panel), 4-site 
with standard patching (second from left), 4-site with alternative patching (4*), (central panel), 8-site (second from right), 
and 16-site (rightmost panel). Momentum-space patches indicated by shaded regions; electron self-energy is independent of 
momentum within a patch but may vary from patch to patch. Dots (red) represent the K points in reciprocal space associated 
to the patches in the DCA construction (see text). Thin lines : Fermi surfaces for the noninteracting system with t' — — 0.15t 
for half filling and hole dopings of 10%, 20%, and 30%. All clusters have an inner patch around (0, 0) (yellow) and an outer 
patch around (n, 7r) (green). Clusters with four or more sites also have an antinodal patch at (n,0) and symmetry-related 
points (blue), clusters with eight or more sites have a nodal patch ((7r/2, 7r/2), red). The 16-site cluster has two additional 
independent momentum sectors, around (7r/2, 0) (orange) and around (37r/2, 7r/2) (cyan). All clusters have the full point-group 
symmetry of the lattice. 



occur, the features of the cluster which cause them are 
determined. Our results define the current frontier of the 
field, given present computational capabilities, and call 
for a new generation of theoretical developments aiming 
at improving momentum-space resolution. 

While the various aspects of the doping-dependent 
phase diagram of the two-dimensional Hubbard model 
have been noted in various ways in the cluster dynami- 
cal mean-field literature, the generality of the results and 
their robustness to choice of cluster have not been previ- 
ously appreciated. The comparison of results for different 
sized clusters clearly demonstrates that the essentials of 
the carrier concentration dependence of physical proper- 
ties of a doped Mott insulator are as sketched in Fig. [T] 
Far from the insulating state, the properties are those of 
a moderately correlated Fermi liquid. Moreover, the mo- 
mentum dependence of the renormalizations is very weak: 
the properties are described well by single-site dynamical 
mean-field theory, as previously noted, e.g., in Refs. [32j 
and |33| . We refer to this regime as the isotropic Fermi 
liquid. [Note that "isotropic" here means isotropic scat- 
tering properties (self-energy) along the Fermi surface, 
but the Fermi surface is not circular.] As the doping is 
decreased towards the n = 1 insulating state the sys- 
tem enters an intermediate doping regime where the low- 
temperature behavior is still described by Fermi-liquid 
theory, but the Fermi liquid is characterized by a strong 
momentum dependence of the self-energies with the self- 
energies being largest near the zone corner (0,7r)/(7r,0) 
points and smallest near the zone diagonal (±7r/2, ±7r/2) 
regions of momentum-space. We refer to this as the 
regime of momentum- space differentiation. The change 
between the isotropic and momentum-space differenti- 
ated Fermi-liquid regimes is not characterized by any or- 
der parameter and we believe it to be a crossover, not a 
transition, but the doping at which the change occurs is 
surprisingly sharply defined and is indicated by dashed 
lines in Fig. [T] 



As the doping is decreased yet further, a non-Fermi- 
liquid regime appears on the hole-doping side but not on 
the electron-doping side (for the moderate particle-hole 
asymmetry considered here). In the non- Fermi-liquid 
regime, regions of momentum-space near (0, ir)/(ir, 0) ac- 
quire an interaction-induced gap, while the zone diagonal 
regions of momentum-space remain gaplcss. We refer to 
this regime as the sector-selective regime. The boundary 
between the regime of momentum-space differentiation 
and the sector-selective regime is indicated by a light 
solid line in Fig. [T] Finally at doping n = 1 the system 
is in the Mott insulating phase. 

The remainder of the paper is organized as follows. In 
Sec. |n] we summarize the general features of the doping- 
driven Mott transition, define the model to be studied 
and the questions to be considered and outline the theo- 
retical approach. In Sec. lIIII we demonstrate the existence 
of different doping regimes and how they appear in the 
different cluster calculations. Section lTVl cxplores in more 
detail the intermediate "momentum-space differentiated" 
doping regime, studies the momentum-selective regime 
and aspects associated with the pseudogap. Section IVl 
then considers the sector-selective regime. In Sec. lVII wc 
summarize our insights into the behavior of smaller size 
clusters. Finally, Sec. I VIII is a summary and conclusion, 
also pointing out directions for future work. 



II. MODEL AND METHOD 

In conventional electronic-structure theory, band insu- 
lators are periodic crystals in which all electronic bands 
are either filled or empty. A necessary condition for band- 
insulating behavior is that the number of electrons per 
unit cell is even. For the purpose of this paper we define 
a correlation-induced or "Mott" insulator as a periodic 
crystal which has no broken symmetry but is insulating, 
even though the number of electrons per unit cell is not 



an even number. The insulating behavior must arise from 
electronic correlations which are beyond the scope of con- 
ventional band theory. The basic theoretical model used 
to study the correlation-induced transition is the one- 
orbital Hubbard model, consisting of electrons hopping 
among sites i of a lattice and subject to an interaction 
U which disfavors configurations with two electrons on a 
given lattice site. The model is most conveniently written 
in a mixed position-spacc/momcntum-spacc representa- 
tion as 



H 



E 



£kC k<y c kcr 



^E 



n^riii. 



(1) 



Here c ka is the Fourier transform of the operator c ia 
which creates an electron of spin a on site i and hi a = 
c ia Ci a is the number operator for electrons of spin a on 
site i. We specialize to a two-dimensional square lattice 
with hopping parameters such that the electron disper- 



se 

e k = — 2£ (cos k x + cos k y ) — At' cos k x cos k y 



(2) 



We choose energy units such that t = 1 and special- 
ize to t' = —0.154. A nonzero t' introduces a particle- 
hole asymmetry into the problem. The relative magni- 
tude and sign of t' arc chosen to be representative of the 
band structure of high temperature superconductors^ 
we note that for very much larger amplitudes of t' 
(\t'\ ~ 0-3) the behavior, especially on the electron- 
doped side, differs slightly from what is discussed here4^ 
We take U = It; this value is slightly smaller than 
the U ~ 9t believed to be relevant to high-temperature 
superconductors^ but is large enough that for all of the 
clusters we study a Mott insulating phase exists while be- 
ing small enough to enable us to perform computations 
on large clusters with the resources available to us. 

To solve the model we use the "dynamical clus- 
ter" (DCA) version of cluster dynamical mean-field 
theoryi 23 i 27 In this approach the Brillouin zone is par- 
titioned into N equal-area tiles which we may label 
by N central momenta K{. The electron self-energy 
E(fc, u>) is approximated by a piccewise constant function: 
S(fc, u>) = Ejq (ui) if k is contained within the tile labeled 
by Ki. The functions T,k i (cj) are obtained from the solu- 
tion of an auxiliary iV-site quantum impurity model with 
on-site interaction U and hopping and bath-coupling 
parameters determined by a self-consistency condition 
which may be written 



G K (iu n ) 



K 



d 2 k 



1 



iu n + n -£fc - T. K (iuj n )' 



(3) 



Here Gk and Y>k are Green's functions and self energies 
of the impurity model and J denotes the normalized 
integral over the momenta in the sector labeled by mo- 
mentum K. For details of the self-consistency procedure 
see Ref. M- 

DCA schemes are distinguished by the number and 
arrangement of the tiles. Note that while (as far as is 



known) the method requires all of the tiles to be non 
overlapping and to have a volume equal to 1/N of the 
Brillouin zone both the shape of the tiles and their ar- 
rangement may be freely chosen. The freedom to adjust 
the tile shape was introduced in recent 2-site dynamical 
mean-field studies i 37 ! 38 We restrict attention to tilings 
which respect the point symmetry of the lattice but con- 
sider different tiling shapes. We use the 2-site patching 
introduced in Ref [33 and introduce an alternative patch- 
ing, which we refer to as 4*, for the 4-site cluster. The 
specific clusters we consider are shown in the five panels 

of Fig. m 

Also shown as light lines in Fig. [2] are the Fermi sur- 
faces of the noninteracting model for carrier concentra- 
tions corresponding to half filling and 10%, 20%, and 
30% hole doping. The Fermi surface lines show a de- 
ficiency of the standard 4-site cluster: for holc-dopings 
near to half filling the Fermi surface is almost entirely 
contained in the (0, n)/(n, 0) sectors; thus this cluster 
has a difficult time capturing momentum-space differ- 
entiation along the Fermi surface. The alternative 4* 
patching shown in the middle panel of Fig. [5] offers the 
possibility of capturing some of the zone-diagonal region 
of the Fermi surface within a different sector. 

The essential computational task is the solution of the 
quantum impurity model. To accomplish this we use 
continuous-time quantum Monte Carlo methods in the 
auxiliary field (CT-AUX) formulation^. This method 
is an imaginary time method which yields the particle 
densities in each sector, along with the sector Green's 
function Gk and the sector self-energy T,k- From Gk 
we obtain the sector occupancy tik via 



n K 



G k (t ->. 0- 



= \iZ G ^ 



(4) 



Note that the sector Green's functions are normalized 
in such a way that for tik = 2 all k states in a sector 
are occupied by two electrons. The total density is n = 
2(J2k n K)/N (the two is for spin). 

In our analysis we work for the most part with sector 
quantities Gk, ^k, and tik- We prefer to avoid the 
"periodization" or interpolation schemes which attempt 
to reconstruct continuous functions of momentum from 
the coarse-grained quantities which are the direct output 
of the calculation. 

Important quantities for the following discussion are 
the parameters fi* k , Z k , and £"(&, 0) = H"(k,u = 0). 
These are defined generally for a Fermi liquid in terms 
of the low-frequency limit of the real (E ) and imaginary 
(S ) parts of the retarded electron self-energy E(fc, w) as 



E(fc, u) « nl - fi + iE"(k, 0) + w(l - Z k l ) + .. 



Zl 



l-d u V(k,u)l 



It will also be useful to consider 

r fc = z fe |s"(fc,o)| 



(5a) 
(5b) 
(5c) 

(6) 



In the DCA approximation we use here these become 
piecewise constant functions of momentum; we denote 
the value appropriate to sector K by suppressing the mo- 
mentum argument and adding a subscript K. 

In the Fermi-liquid regime, these parameters express 
important aspects of electronic physics. For completeness 
we briefly recall their meaning here. At low frequencies 
ui — >• 0, the spectral function A(k, w) = --ImG(fc,w) 
becomes 



fromi£ 



A(Je,w) 



1 



Z k T k 



7T {uJ + z k ^l-e k )) 2 +rl 



(7) 



fifc determines the location of the renormalized Fermi sur- 
face (which is the locus of points fci? for which e kp = 
Hj, ). Thus a momentum dependence of [i* k signals a 
change in shape of the Fermi surface and more gener- 
ally a shift in the mean energy of one momentum sector 
relative to the others. T k is the width of the quasipar- 
ticle peak. For a fixed k on the Fermi surface, A(kp,uj) 
is peaked at w = and the width in frequency is set 
by T k . A necessary condition for Fermi- liquid behavior 
is that T k be small, in which case A(k,u) is character- 
ized by a reasonably well-defined quasiparticle peak with 
frequency width given by T k and area given by the quasi- 
particle weight Z k . The criterion T k < itT is the math- 
ematical expression of the condition that the width of a 
thermally excited quasiparticle is less than its energy. In 
a Fermi liquid E"(fc,0) - T 2 as T -> so T is para- 
metrically less than T. In this paper, however, since we 
cannot reach very low temperatures, it will be useful to 
relax this definition and consider as a quasi-Fermi liquid 
any system where (i) the Luttinger theorem is reason- 
ably well obeyed, (ii) T k decreases as T decreases at all 
points along the renormalized Fermi surface, and (iii) at 
all points along the Fermi surface T k < ttT. 

Defining the bare velocity v k — de k /dk, the disper- 
sion away from the Fermi surface is determined by the 
renormalized velocity 



v k 



Z k 



de k 
dk 



dX'(k,w = 0) 
dk 



Z K v k , 



(8) 



where in the second equality we used the fact that in the 
DCA approximation used here the self-energy is a piece- 
wise constant function of momentum so that the quasi- 
particle weight Zk coincides with the velocity renormal- 
ization. 

Note that dc transport properties are controlled by the 
mean free path 



Z fc = E (fc jW )/ dts k + d^(k,u = 0) 



X k (lu = 0) 



''/.■ 



(9) 

i.e., by T,"(K, 0), not by Tk- We have used the piecewise 

continuity of the DCA self-energy in the second equality. 

An estimate for the sector spectral function Ak(oj) in a 

range of order T around the Fermi level may be obtained 



('f 



dio Ak(ui) 



2T cosh 



(10) 



2T 



because the cosh - term implies that this integral is 
sharply peaked around u> = 0. 

The quantities [i* K , Zk and E^(u; = 0) are defined 
in terms of the real-axis self-energy while the impurity 
models are solved on the imaginary frequency axis. Be- 
cause the studies in this paper require the survey of a 
wide range of clusters and dopings we have not been able 
to obtain data of sufficient quality to permit reliable an- 
alytical continuation of all of our cluster data. Different 
size clusters impose different computational burdens and 
(with the computational resources we have) lead to data 
with differing statistical errors, introducing further un- 
certainty in the comparison of continuations based on 
data from differently sized clusters. We therefore esti- 
mate \x* K , Zk and the spectral function at zero energy 
Ak{oj = 0) directly from the Matsubara-axis self ener- 
gies. n K is obtained by extrapolating the real part of the 
Matsubara-axis self-energy to zero. E^-(w = 0) may be 
obtained as the extrapolation of the imaginary part of the 
Matsubara self-energy to zero frequency. The Kramers- 
Kronig relation connecting E and E implies that if E" 
is small at low frequencies (as is the case in a Fermi liq- 
uid) Zk can be obtained from the extrapolation to zero 
frequency of the derivative of the Matsubara-axis E(iw„). 
We therefore estimate n* K , Zk, and Y,"{K, 0) by fitting 
the three lowest calculated Matsubara frequencies to a 
quadratic form, which we use to extrapolate the value 
and first derivative to zero Matsubara frequency. 

This "poor man's analytical continuation" becomes 
less reliable at higher temperatures because the spac- 
ing between Matsubara points is 2itT and therefore dis- 
cretization errors become larger for higher T. At all 
temperatures we study we find that of the three quan- 
tities fi K nas t ne least uncertainties (in the Fermi-liquid 
regime) because the real part is continuous across the 
real-frequency axis and has a relatively weak temperature 
dependence. H"(K, 0) is less accurately determined than 
fi* K because it has a marked temperature dependence and 
becomes small at low T, and Zk is least accurately de- 
termined because it is a derivative. 

For momenta removed from the Fermi surface the in- 
teresting physics (in a Fermi liquid) requires knowledge of 
the self-energy at nonzero frequency, in particular at the 
frequency u> qp which solves the "quasiparticle equation" 
u; gp + E'(fc, uj qp ) — e k — 0. Our poor man's analytical con- 
tinuation procedure does not provide much information 
about the real-frequency self-energy at non-zero frequen- 
cies, and the quantities n K , Zk, and Ti"(K, 0) do not 
have a clear meaning for states with e k ^ /i* . We there- 
fore focus in the following on momentum sectors which 
contain pieces of the Fermi surface and on regimes where 
reasonably long-lived quasiparticles exist. 

Further comparison of our methods to analytic contin- 
uation is given in appendix |A"1 
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TABLE I: Characteristic densities for onset of momentum dif- 
ferentiation on the hole-doped (n^ iS ) and electron-doped sides 
(n^iff), estimated gap of the Mott insulator (A 9 ), location of 
the sector-selective transition on the hole-doped (n§ ST ) side 
(there is no SST transition on the electron-doped side for 
the parameters chosen in this paper), and gap of the sector- 
selective regime (Asst). 



FIG. 3: Main panel: total electron density n as a function of 
chemical potential jj, for clusters considered in this paper at 
inverse temperature fit = 20 for all clusters except 16, where 
fit = 7.5 is shown. Inset: expanded view of small chemical- 
potential region, highlighting region of Mott gap where n — 1 
independent of \i. 



III. DIFFERENT DOPING REGIMES 

In this section we characterize the different doping 
regimes discussed in the introduction in terms of the de- 
pendence of the total electronic density n on the chemi- 
cal potential fj, and in terms of the partial occupancy Hk 
of the patch associated with momentum K, defined by 
Eq. ®. 

The total electron density n is plotted in Fig. [3] as 
a function of the chemical potential [i, for all clusters 
studied in this paper. For comparison, the single-site 
dynamical mean field (DMFT) result (i.e., the cluster 
with N = 1) is also displayed. Data are shown at the 
inverse temperature fi = 20/i for all clusters except 16 
sites, where we used (3 = 7.5/i. We first observe that at 
half- filling (n = 1) all N > 1 cluster calculations yield 
an insulating state, revealed by the plateau in the n(/i) 
curve. By contrast, single-site DMFT yields a metallic 
state [no plateau in n(//)], because the U/t = 7 stud- 
ied here is smaller than the single-site DMFT value of 
the critical interaction strength U « 12£. The difference 
occurs because cluster calculations take spatial correla- 
tions into account, and these stabilize the insulating state 
at smaller U values than are needed in the single-site 
calculation^^^SiSi 

The range in chemical potential over which the n(/j,) 
curve is flat can be used to define an estimate A g for 
the insulating gap as A g = ji(n = 1 + ) — ji(n = 1~). Of 
course thermal effects mean that the n([i) curve is not 
precisely flat. We adopt the criterion that the chemical 
potential is within the gap if 0.99 < n < 1.01. Table U 
presents results obtained for the different clusters, using 
this criterion, at inverse temperature f3t = 20 (fit = 7.5 
for 16 site). Except for the conventionally patched 4- 
site cluster, all of the clusters studied give remarkably 
consistent estimates for A„. We will explain below why 



the conventionally patched 4-site cluster is different. 

We next turn to the doping dependence. It has pre- 
viously been shown from small cluster studies that at 
large enough doping, cluster corrections to the single-site 
DMFT results are smalli^i^ Figure [3] shows that both on 
the hole-doped and electron-doped sides the differences 
between the results for different clusters become small 
so that the n(^i) curve is well described by the single- 
site DMFT result, consistent with previous findings. Be- 
low a critical chemical potential or critical density the 
curves separate, indicating dependence on cluster geom- 
etry and therefore momentum differentiation. We define 
critical carrier concentrations n'^ iB (on the hole-doping 
side) and n^g (on the electron-doping side) at which a 
significant momentum dependence appears in the elec- 
tron self-energy. The ndis may be estimated from the 
densities at which the n(/x) curves begin to separate, 
but as we shall see a better estimate may be obtained 
from analysis of the sector dependence of the renormal- 
ized chemical potential \i* K as discussed in Sec. IIV1 It is 
this latter analysis which is used to obtain the densities 
shown in Table [U 

Figures 2] and [5] display the occupancies tik of the dif- 
ferent sectors, respectively, as a function of the chemical 
potential fj, and the total density n. For the clusters with 
N < 16 these figures display data for two temperatures: 
[3t = 20 (symbols and heavy lines) and fit — 7.5 (light 
lines), showing that the temperature dependence is weak. 

Figure 0] shows that for all clusters except the 2 site, 
there is a range of chemical potentials where the density 
in the (0,7r) sector remains locked at n^^) — 1/2 cor- 
responding to a half-filled patch. The chemical-potential 
range over which the density in the (0, 7r) sector is flat 
defines a gap, Asst, which is also shown in Tableland is 
of the same order for all clusters. The presence of a gap 
is confirmed by the density of states estimates discussed 
in Sec. El 

The regime of chemical potentials where the (0, 7r) sec- 
tor is incompressible while other sectors can be doped is a 
signature of the sector-selective regime discussed in pre- 
vious work i 40 ' 41 We see that the phenomenon is robust, 
occurring in 4*-, 8-, and 16- site clusters. It is also inter- 
esting to note that while (for the parameters we studied) 
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FIG. 4: Partial occupancies uk in the different patches of the Brillouin zone, as a function of chemical potential /i, for different 
clusters. Inverse temperature is j3t = 20, thin lines and 16-site cluster: fit = 7.5. The arrows indicate the points where 
the sector densities deviate from the plateau values; the chemical-potential difference between the arrows defines the range of 
chemical potential corresponding to the Mott insulator (A 9 ) and to the sector-selective regime (Asst)- The Mott insulating 
density n — 1/2 is shown as the dotted line. 




FIG. 5: Partial occupancies uk in the different patches of the Brillouin zone, as a function of total electron density n for the 
different clusters. Inverse temperature fit = 20, thin lines and 16-site cluster: fit = 7.5. 



the conventionally patched 4-site cluster does not have a 
sector-selective regime, the change in momentum-space 
patching involved in going from 4 to 4* gives the clus- 
ter the possibility of distinguishing nodal and antinodal 
excitations, and the transition reappears. This, along 
with the similar behavior of the 8- and 16-site clusters, is 
strong evidence that the sector-selective transition (SST) 
is a generic phenomenon that appears wherever the mo- 
mentum patching allows it. It reveals strong momentum 



differentiation and the formation of an antinodal pseudo- 
gap at low hole-doping. 

A different perspective on the (0, 7r) sector gap is pro- 
vided by the plots as a function of the total density shown 
in Fig. [5] We see that in the 4*-, 8- and 16-site clusters 
the partial density in the (0,7r) sector remains locked at 
n (o.7r) = 1/2 for a range of densities while the density 
in other sectors changes. This behavior, that for some 
range of doping the antinodal patch remains incomprcss- 



ible while the other patches accommodate the dopants, 
is observed only on the hole-doped side (for the parame- 
ters U/t, t'/t studied here) and defines the sector-selective 



regime. The boundary, 



l SST' 



of the density range 



over which the (0, 7r) sector is incompressible is a true 
T = phase transition of the self-consistent DCA equa- 
tions associated with a specific cluster. For concreteness 
we define rigg T as the value of the total density n, where 
\ n (o 7r) — 0.5 1 = 0.005 and we report these values in Ta- 
ble ID 

Unlike the 4*-, 8- and 16-site clusters, the conven- 
tionally patched 4-site cluster does not display a sector- 
selective regime. However, we see from Table U that the 
insulating gap for the conventionally patched 4-site clus- 
ter has a magnitude very similar to the (0, 7r) sector gap 
found in the 4*-, 8- and 16-site clusters. This makes 
the origin of the difficulties of the 4-site cluster clear: in 
this cluster the (0,7r) patch must fill two roles. First, 
this patch contains the (0,7r) point so that its behavior 
must represent the physics of the pseudogap. Second, 
the patch contains essentially all of the Fermi surface, 
and therefore must represent the behavior of the gapless 
nodal quasiparticles. The modified 4* patching avoids 
this problem because the (0, 0) sector contains the nodal 
portion of the Fermi surface so that in this cluster the 
(0, 7r) sector does not have two roles. We also note that 
the critical density for the onset of the sector-selective 
regime is about 0.92 for the 8- and 16-site clusters, but 
closer to n = 1 for the 4* cluster. We believe the differ- 
ence arises because, while the 4* cluster does allow some 
nodal/ antinodal differentiation, the portion of noninter- 
acting Fermi surface which is outside of the (0, n) sector 
remains comparatively small. 

The 2-site cluster is also different: neither of the two 
momentum sectors exhibits a plateau outside of the insu- 
lating regime although previous studies of the two-patch 
cluster using a larger t' (Refs. [37] and 1381 ) have revealed 
a sector-selective regime for a range of hole-dopings. We 
shall see in Sec. [V] that the difference arises from a dif- 
ference in the nature of the sector-selective transition in 
this cluster. 



IV. MOMENTUM-SPACE DIFFERENTIATION 

In this section we analyze in more detail the 
momentum-space differentiated regime introduced in the 
previous section. We begin our discussion with \x* K , 
Eq. (|5b[) . shown as a function of total density n in Fig. |6] 
For sufficiently high levels of electron or hole doping, n* K 
becomes almost independent of momentum sector within 
a cluster and indeed takes the same value independent of 
cluster. As the doping is decreased toward half filling the 
traces separate, indicating that the self-energy begins to 
depend on momentum. The onset of momentum-space 
differentiation does not (to our knowledge) correspond to 
a phase transition. Indeed at any carrier concentration 
the self-energy has some momentum dependence: per- 
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FIG. 6: (Color online) Renormalized chemical potential [i* K = 
/i — T,'(K, uj — 0) for different momentum sectors, plotted vs. 
particle number n in (from top to bottom panels) 2-, 4-, 4*-, 
8- and 16- site clusters. Dotted line (green): location of lower 
band edge for sector (n, n) in the 2-site cluster. 



haps weak but, in general, nonzero so that any definition 
of the doping at which the variation with momentum 
becomes significant involves some arbitrariness. In Ta- 
ble HI we estimate this density n^ iff (on the hole-doping 
side) and n^ s (on the electron-doping side) as the density 
where m.&yiK,K' iMif ~ Mlcl > 0-2. 

A momentum dependence of fj,* implies a change in 
shape of the Fermi surface. The two panels of Fig. [7]show 
the Fermi surface computed for the 8-site cluster in one 
Brillouin zone quadrant at two dopings, one (left panel) 
in the high-doping isotropic Fermi-liquid regime and one 
(right panel) in the intermediate doping momentum- 
space differentiated regime. In the isotropic Fermi-liquid 
regime the very weak k dependence of £ means that 
the Fermi surface of the interacting model coincides with 
the noninteracting Fermi surface (at the same density). 




FIG. 7: Plot of the upper quadrant of the Brillouin zone show- 
ing the Fermi surface as a black line. The dashed line shows 
the noninteracting Fermi surface with the same total density. 
The sectors of the 8-site cluster are also displayed (see Fig. 01 ■ 
Left panel: n » 0.69 (in isotropic Fermi-liquid regime); right 
panel n ~ 0.82 (in momentum-space-differentiated regime). 
All Fermi surfaces are consistent with the Luttinger theorem. 



However, in the momentum-space differentiated regime 
the shape changes. The piecewise constant nature of 
the DCA self-energy means that the DCA Fermi surface 
is also a piecewise constant approximation to the true, 
renormalized Fermi surface. We see that the interaction 
has the effect of enhancing the Fermi-surface curvature 
so that the renormalized Fermi surface corresponds, in 
effect, to a quasiparticle band structure with a bigger t' 
than the bare theory, consistent with the suggestions of 
Refs.[Hand[H 

An important result of Fermi-liquid theory is the "Lut- 
tinger theorem,"— which states that the volume of k 
space contained within the Fermi surface, defined by 
Ek + Re[E(fc, u = 0)] = yU, is equal to the density of par- 
ticles per unit cell (mod2). We have verified that both 
Fermi surfaces are consistent with the Luttinger theorem 
to within 1%. Corresponding plots made for dopings out- 
side the Fermi-liquid regime are not consistent with the 
Luttinger theorem because the "Fermi surface," although 
mathematically defined, is not physically meaningful. We 
also note that the proof of the Luttinger theorem relies in 
an essential way on an integration over the whole zone; 
it is not obeyed sector by sector, meaning that the tik 
defined from the sector Green's function need not cor- 
respond to the sector area contained within the Fermi 
surface. 

It is worth emphasizing that the isotropic Fermi-liquid 
regime, in which single-site DMFT is barely corrected by 
cluster expansions, is not a trivial weak-coupling regime. 
This may be seen from Fig. [5J which shows the full fre- 
quency dependence of the imaginary parts of all com- 
ponents of the self-energy of the 8-site cluster for dop- 
ings n = 0.72 (in the isotropic Fermi-liquid regime) 
and n = 0.9 (momentum-differentiated regime). In the 
isotropic regime the relative differences in self-energy be- 
tween different sectors arc small whereas they are large 
in the differentiated regime. 

We next turn to the scattering rate for low-energy ex- 
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FIG. 8: Imaginary part of Matsubara self-energies of the 8- 
site cluster for two densities at inverse temperature fit — 20, 
for various sectors. 



FIG. 9: | E" (if, 0)| for different momentum sectors plotted vs 
particle number n in (from top to bottom panels) 2-, 4-, 4*-, 
8-, and 16- site clusters. 



citations, which we estimate from T*"(K, 0), the extrapo- 
lation to iuj n = of the imaginary part of the Matsubara 
self-energy [Eq. ([5a])]. Figure [§J presents results obtained 
at the relatively high temperature fit = 7.5 where we 
have data for all clusters and dopings. We see that at 
high dopings the scattering rates are small and (within 
our numerical accuracy) independent of cluster size and 
sector, reflecting the reasonably isotropic nature of the 
high-doping regimes. As doping is decreased the curves 
clearly exhibit an onset of momentum-space differentia- 
tion at a doping ~ 0.25 consistent within errors with the 
doping estimated from fi* and again only weakly depen- 
dent on cluster size. 

Momentum-space differentiation is marked by a strong 
relative increase with decreasing doping of E"((0, tt),0) 
and a less strong increase in £"((7r/2,7r/2), 0) (for the 
clusters which provide access to this sector). These are 
the two sectors which contain the Fermi surface and for 
which T,"(K, 0) has meaning as a scattering rate. Com- 
parison of electron and hole dopings shows that while 
momentum-space differentiation sets in at about the 
same absolute value of doping in the two cases, the degree 
of differentiation between sectors (0,7r) and (tt/2,tt/2) is 
greater on the hole-doped side than on the electron-doped 
side. 

We turn now to a more detailed examination of results 
from the 8-site cluster, which is large enough allow a di- 
rect comparison of the nodal and antinodal regions of 
the Fermi surface, but is small enough to allow detailed 
computations down to relatively low temperatures. The 
two panels of Fig. [TU] show H"(K,0) and the quasipar- 
ticle weight/velocity rcnormalization Zk for the nodal 
K = (it/2, ir/2) and antinodal K = (%, 0) sectors as a 
function of doping at a relatively low and a relatively 
high temperature. 

Comparison of the two panels of Fig. [TU] shows that 
the momentum-space differentiation is marked primar- 



10 



ily by a variation in scattering rate. As doping is re- 
duced, the Fermi-surface scattering rates increase rapidly 
and a marked difference between the two Fermi surface 
sectors develops with the antinodal sector K = (ir, 0) 
characterized by a much more rapidly growing scatter- 
ing rate. Further, the scattering rates exhibit a pro- 
nounced particle-hole asymmetry. However, while the 
inverse mass enhancement/velocity rcnormalization Zk 
decreases as doping is decreased, the variation with dop- 
ing is much less dramatic and, interestingly, there is very 
little particle-hole asymmetry or difference between the 
two momentum sectors. We also note that the nodal 
quasiparticle residue Z( 7r / 2 ,n/2) appears to extrapolate 
to a non-zero value at n = 1. (A different result was 
found using self-energy interpolations in superconducting 
state cellular dynamical mean field (CDMFT) calcula- 
tions on 4-site clusters . 58 i 59 ) This is inconsistent with the 
Brinkman-Rice theory but qualitatively consistent with 
data on high-T c materials, where photoemission mea- 
surements indicate a zone-diagonal quasiparticle velocity 
which is only weakly doping dependent^ (Very recent 
measurements indicate that if the velocity is measured 
on very low scales, below the resolution of the numerics 
in this paper or of previous photoemission data a stronger 
doping dependence of the velocity is found^) 

Figure [TT] presents the temperature dependence of the 
nodal and antinodal scattering rates obtained for the 8- 
site cluster for selected densities. To highlight the tem- 
perature dependence we plot £ jT . While the temper- 
ature range accessible to us is too limited to establish 
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FIG. 10: Estimators of quasiparticle weight /velocity renor- 
malization factor Z (upper panel) and scattering rate 
|£" (K, 0) | obtained from extrapolation of Matsubara-axis self 
energies for the K = (n, 0) and K = (tt/2, n/2) sectors of 8- 
site cluster at temperatures T = 4/7.5 and t/20, along with 
estimates obtained by maximum entropy analytical contin- 
uation (a.c). Vertical dashed lines denote the boundaries 
separating the regime consistent with Fermi-liquid theory 
(/3r_R- < 7r) from the regime that is clearly not Fermi liquid. 
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FIG. 11: (Color online) Temperature dependence of 
| E" (if, 0)| for (-7T, 0) sector (dashed lines, black) and 
(7r/2, 7r/2) sector (solid lines, red) of 8-site cluster at dop- 
ing 0.7 (open symbols), 0.75 (light symbols), and 0.8 (filled 
symbols). 



any specific form of temperature dependence it is clear 
that at the higher doping isotropic Fermi-liquid regime 
(n = 0.70), the scattering rates drop faster than linearly 
at low T while at the lower doping (n = 0.80; momentum- 
differentiation regime) the two sectors have different tem- 
perature dependence at low temperature with the nodal 
sector scattering rate vanishing more rapidly than T 
at low T and the antinodal rate vanishing less rapidly. 
At the intermediate doping n = 0.75 on the bound- 
ary between the two regimes the behavior is intermedi- 
ate. These features arc in qualitative agreement with 
the momentum-space variation in the electronic mean 
free path inferred from angular-dependent magnetore- 
sistance experiments*^ For highly overdoped cupratcs 
these experiments reveal a scattering rate which is rea- 
sonably isotropic around the Fermi surface and exhibits 
a relatively conventional temperature dependence. Be- 
low a critical doping a momentum-space differentiation 
appears with the antinodal scattering rate being larger 
and exhibiting a weaker temperature dependence. 



V. SECTOR SELECTIVE REGIME 

In this section we discuss in more detail the sector- 
selective regime, which has been characterized in Sec. IIIII 
by the existence of a plateau in the uk(^) curve for the 
sectors containing (0, 7r) (and symmetry-related points) 
and the absence of such a plateau in the sectors con- 
taining (±7r/2, ±7t/2). We first observe that, as seen, 
for example, in Fig. O the approach to the sector- 
selective phase in clusters 4*, 8, and 16 is marked by a 
rapid increase in the absolute value of E"((0, ir), 0) and 
£"((7r/2, 7r/2),0), with E"((0,7r),0) being much larger. 
In the (0, it) sector the increase is associated with the 
formation of a pole in the sector self-energy. This pole 
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is responsible for the sector-selective insulating behav- 
ior, as discussed in detail in Refs. |40h421 However, we 
also see that as doping is decreased the (tt/2, 7r/2) sector 
self-energy increases, and in fact becomes large enough 
that (at least at the temperatures accessible to us in this 
study) this sector is clearly in a non-Fermi-liquid regime 
(although as seen it is compressible). It is an interesting, 
but so far unresolved, question whether the nodal sector 
is intrinsically non- Fermi liquid in this doping regime or 
whether the Fermi temperature is simply lower than the 
lowest temperatures we can access. 

We next show that, as expected on physical grounds, 
the plateau in sector density tik coexists with a gap in 
the sector density of states. Figure [T2l shows (3\Gk(P/2)\ 
[Eq. [10] which, as T — > 0, converges to the sector density 
of states at the Fermi level, Ak(u = 0). We first fo- 
cus on the lower left panel, showing results for the 8-site 
cluster at several temperatures. We see, as previously 
shown , 40 ! 41 that the (0, 0) and (tt, tt) sectors have essen- 
tially no low-energy density of states. This is expected 
because the Fermi surface does not pass through these 
points at any of the densities displayed. The (tt/2, ir/2) 
sector shows a density of states which at higher doping, 
in fact, takes the Fermi-liquid value, but in the sector- 
selective regime is suppressed from the Fermi-liquid value 
[consistent with the large value in this regime of the scat- 
tering rate in the (tt/2, tt/2) sector] but at lower temper- 
atures is weakly temperature dependent. The weak tem- 
perature dependence suggests, but with our limited tem- 
perature range docs not prove, that the sector is in a non- 
Fermi-liquid regime. By contrast the sector (0, tt) results 
display a strong temperature dependence: for carrier con- 
centrations above ~ 0.9 the sector density of states is a 
strongly decreasing function as temperature is decreased, 
indicating that a gap is forming in this sector. 

Now consider the other panels in the plot, beginning 
with the 16-site cluster. The temperature range accessi- 
ble to us is more limited, because at high doping our abil- 
ity to calculate at low temperature is hampered by a sign 
problem. However, within the available resolution we see 
that the 16-site cluster behaves in the same manner as 
the 8-site cluster. This is in contrast to the convention- 
ally patched 4-site cluster, where we see that the density 
of states in the (0, tt) sector, while suppressed and weakly 
temperature dependent, only vanishes at n — 1, similar 
to the behavior of the (ir/2, tt/2) sector in the 8- and 16- 
site clusters. The (0, 0) sector is simply gapped, similar 
to the (0, 0) sectors in the larger clusters. The alterna- 
tive patching of the 4* cluster remedies this behavior by 
including some of the nonintcracting Fermi surface in the 
(0, 0) sector, and we see that this cluster behaves more 
similarly to the 8- and 16-site clusters. In the 4* cluster 
the sector-selective transition exists but is seen at a lower 
doping than in the larger clusters [corresponding to the 
smaller part of the Fermi surface that (0,0) intersects], 
and the behavior of the (0, 0) sector is intermediate be- 
tween that of a nodal sector in a larger cluster and the 
conventionally patched 4-site cluster. 



VI. SMALL VERSUS LARGE CLUSTERS 

In this section we compare the results obtained on 
small (2- and 4-site) clusters to those obtained on the 
larger 8- and 16-site clusters. Roughly speaking, the sec- 
tor quantities Gk and Ejf correspond to averages of the 
intrinsic quantities G(k,uj) and Y>(k,uj) over appropri- 
ate regions of the Brillouin zone. As the cluster size 
is reduced, the averages extend over wider regions of 
momentum-space. While it is clear that very small clus- 
ters will only be able to capture a very coarse-grained ver- 
sion of momentum-dependent physics, it is nevertheless 
interesting to understand the strengths and weaknesses 
of small cluster studies, both because it is probably not 
feasible to use larger clusters to study multiorbital mod- 
els and because even for the one orbital Hubbard model 
much more accurate data can be obtained, down to much 
lower temperatures, on small clusters than on large clus- 
ters. 

We begin by reconsidering the cluster size dependence 
of the scattering rates shown in Fig. [9l starting from the 
16-site cluster down to the smallest clusters. The self 
energies obtained on the 16- and 8-site clusters for the 
K = (0, 0), (n, 7r), (0, tt), and (tt/2, tt/2) patches are very 
similar. In both of these clusters the different patches 
correspond to regions of the zone with different physi- 
cal properties. When the momentum resolution starts 
to be reduced, the patches become bigger and the clus- 
ter quantities involve averages over regimes with differ- 
ent physical properties and therefore yield results which 
are in some sense intermediate. For example, the (ir, 0) 
sector of the 4-site cluster is seen to have an intermedi- 
ate behavior between the (tt, 0) and (tt/2, tt/2) sectors of 
the 8- or 16-site clusters while the (tt, it) sector of the 
4* cluster has behavior somewhat closer to that of the 
(tt, 0) sector of the larger cluster. In general, it appears 
that if this basic fact is borne in mind as results are in- 
terpreted, the smaller clusters do remarkably well in cap- 
turing the important qualitative behavior, although the 
clear quantitative differences call into question the valid- 
ity of the various interpolation schemes. The example of 
the 4 and 4* clusters also shows how a physically moti- 
vated choice of patching (or even a comparison between 
different patching choices) can reveal important physics, 
in this case better separating the nodal and antinodal 
sectors. 

Let us comment in more detail on the specific case 
of the 2-site cluster, which was analyzed in great detail 
in Refs. |37| and |38|. In these papers the second neigh- 
bor hopping t' = — 0.3t was used and a doping-driven 
sector-selective transition was found, but of a different 
character than the sector-selective transitions discussed 
here and in Refs. |4CJ| and |4j]. The transition found in 
Refs. [33 and |38J is caused by a renormalization of the 
effective low-energy chemical potential by the interac- 
tions and can be understood as follows. For a given DCA 
patching of the Brillouin zone, a necessary condition for 
the Fermi surface of the interacting system to penetrate 
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FIG. 12: Estimator for the sector density of states, [3\Gk(P/2)\, as a function of density (hole-doped case only), for temperatures 
ranging from fit — 7.5 to f3t = 25, and for the clusters 4, 4*, 8, and 16. Symbol styles denote patches, line styles denote different 
temperatures. A clear crossing point as a function of temperature is visible in sector (n,0), for the clusters 4*, 8, and 16, but 
is absent for the standard 4-site cluster. 



within the momentum- i^ patch of the Brillouin zone is 
that the renormalized chemical potential \x* K satisfies the 
condition \f K g [e£ ; „, e^LJ, where e^L „,„„ are the min- 
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imum/maximum energy spanned by the noninteracting 
band e^ as k is varied within patch K. This relation 
directly follows from Eq. ((3]) relating the coarse-grained 
Green's function Gk to the self-energy £#-. 

For very high hole doping (low carrier concentration) 
the Fermi surface is contained entirely within the (0, 0) 
patch of the 2-site cluster. As the doping is decreased, 
the Fermi surface expands and in the noninteracting 
model will penetrate the outer patch at some critical dop- 
ing, which depends on £'. However, in the momentum- 
differentiated regime the sector dependence of the renor- 
malized chemical potential in the 2-site cluster is such 
as to push the renormalized Fermi level away from the 
outer patch [in contrast to the larger clusters, where the 
\i* K is such as to push the Fermi surface in the (0, n) patch 
closer to half filling] . This effect gets larger as half filling 
is approached. For the t' and U studied in Refs. |37| and 
|38| . for some range of dopings the Fermi level was able to 
penetrate into the outer patch, but as doping was further 
reduced the increase in [i* K for the outer patch was rapid 
enough to push the Fermi level out again, resulting in 
a sector-selective transition of the band-insulator type. 
For the smaller t studied here the increase in fi* K for the 
outer patch is such as to prevent the Fermi level from 
ever reaching the outer patch. We may therefore regard 



the 2-site cluster as being in a sector-selective regime for 
all hole dopings, where the antinodal states are gapped 
and suppressed compared to the non-interacting system. 
This physics is p articular to the 2-site cluster, but as 
noted in Refs. 1371 andl38l does capture (as best this cluster 
can) the momentum-space differentiation physics which 
is important for cuprates. 



VII. CONCLUSION 

This paper has presented "DCA" cluster dynamical 
mean-field calculations based on 2-, 4-, 8- and 16- site 
clusters of the two-dimensional Hubbard model, with 
hopping parameters appropriate to high-T c materials and 
an interaction which is slightly weaker than in the ac- 
tual compounds, but strong enough to display the es- 
sential physics. The comparison of results obtained 
on different sized clusters shows unambiguously that 
in the two-dimensional Hubbard model with underly- 
ing particlc-holc asymmetry appropriate to high-T c ma- 
terials, the doping-driven metal-insulator transition pro- 
ceeds through several intermediate regimes, but is dif- 
ferent for electron-doped than for hole-doped materials. 
For both hole-doped and electron-doped materials a high 
doping isotropic Fermi-liquid regime gives way as doping 
is reduced to a momentum-spacc-diffcrcntiatcd Fcrmi- 
liquid regime, in which the quasiparticle residue (cquiv- 
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alent to the velocity renormalization in the approxima- 
tion used here) depends only weakly on momentum but 
the electron-scattering rate acquires a strong dependence, 
being small for states near the zone-diagonal and large 
for states near the zone face. The momentum depen- 
dence is less pronounced on the electron-doped side than 
on the hole-doped side. As doping is further reduced 
on the hole-doped but not on the electron-doped side, 
the momentum-spacc-differcntiated Fermi-liquid regime 
gives way to a sector-selective pseudogapped regime, in 
which the states near (0, 7r) are gapped but the states 
near the zone diagonal are not. This sequence of dop- 
ing dependencies, and in particular the nodal/antinodal 
differentiation and the momentum selective pscudogap, 
have for many years been understood to be key fea- 
tures exhibited by the high-T c copper-oxide-based su- 
perconductors. Cluster dynamical mean-field studies 
over the years have uncovered various aspects of this 
behavior. The systematic survey of different dopings 
and cluster sizes presented here has firmly established 
the robustness of this physical picture. Nodal/antinodal 
momentum-space differentiation and momentum selec- 
tivity are robust features that persist independent of clus- 
ter geometry, and there is a remarkably reliable quanti- 
tative agreement across clusters for many of the observ- 
ables presented here. We conclude that this physics is 
a consequence of short-ranged correlations in the two- 
dimensional Hubbard model at intermediate to strong 
coupling. 

Our comparison also allows one to understand the 
strengths and weaknesses of small cluster calculations. 
In particular, we have shown that a simple deformation 
of the standard 4-site cluster designed to better separate 
the nodal and the antinodal region is sufficient to recover 
the sector-selective transition found in 8- and 16-site clus- 
ters. We expect that the comparative information will be 
important for future use of cluster methods in more com- 
plex (multiorbital) systems, where small clusters are the 
only accessible ones. 

This paper has used the "DCA" formulation of clus- 
ter dynamical mean field theory. A different formulation, 
the "Cellular Dynamical Mean Field Theory" (CDMFT) 
has also been introduced^. CDMFT methods have 
not yet been applied to the range of cluster sizes and 
carrier concentrations needed to carry out comparisons 
such as those we have presented here. Recent CDMFT 
papers^— have reported results on momentum space 
differentiation and sector selective transitions providing 
a different (but related) perspective on the issues raised 
in this paper. An extension of the CDMFT results to a 
wider range of cluster sizes and dopings and a systematic 
comparison of CDMFT results to DCA results would be 
very desirable. 

Currently, the two main limitations of cluster meth- 
ods are the limited momentum-space resolution and the 
fcrmionic sign problem. These limit the accessible clus- 
ter sizes and temperatures and currently prevent us from 
reaching a definite conclusion, for example, on the pre- 



cise nature of the low temperature metallic state in the 
nodal region in the sector-selective phase at low tem- 
peratures or on the quantitative evolution of the Fermi 
arcs as a function of doping and temperature. Also while 
the methods give direct access to one-electron proper- 
ties, wide classes of experiments involve "two-particle" 
probes whose analysis requires vertex corrections. New 
theoretical developments are required to overcome these 
limitations. 
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Appendix A: Analytic Continuations 




» pt = 7.5, (0,0) 
• (JtJt) 

1(71,0) 

^ (7C,7l/2) 
(71/2,71/2) 

(71/2,0) 



o 

i. 

s 



-l l 



FIG. 13: lmE(uj n — 0) as a function of density, for all sectors 
of 16-site cluster, at inverse temperature ft = 7.5/t. Filled 
symbols; extrapolation based on three lowest Matsubara fre- 
quencies; open symbols; extrapolation based on two lowest 
Matsubara frequencies. 

In this appendix we compare the results of the "poor 
man's" analytical continuation procedure used in the 
main text to maximum-entropy analytical continuation 
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results obtained using the methods of Ref. [63. We focus 
on the 16-site cluster, where the restriction to high tem- 
peratures forced on us by the limits of our computational 
resources makes the continuations most challenging. The 
solid points in Fig. 1X31 show Y,"(K, 0) for the 16-site clus- 
ter as a function of density, obtained by extrapolation 
from a fit to the three lowest Matsubara points as de- 
scribed in the text. The open points show values obtained 
by extrapolation from the two lowest Matsubara points. 
The strong divergence associated with the (0, tt) sector is 
evident, as is the weaker divergence of the (tt/2, tt/2) sec- 
tor. We would like to interpret these data as indicating 
the presence of a growing scattering rate. However, we 
also see that over a wide doping range the (ir, tt/2) sec- 
tor S"(A', 0) is, in fact, larger than the Z"(K,0) of the 
(tt/2, tt/2) sector, even though the Fermi surface does 
not pass through this sector so that one might expect 
that low-energy processes involving this sector are sup- 
pressed. Of course, for this sector it is not correct to 
interpret T,"(K, 0) as a scattering rate, but the question 
remains of what is the meaning of the result and whether 
it casts doubt on the interpretation of the T,"(K, 0) in the 
sectors that contain the Fermi surface. 
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FIG. 14: Matsubara-axis self-energy for doping x 
solid symbols, lmE(iu n ); open symbols Re£(iu;n). 



0.058; 



To gain further insight we display in Fig Q3]thc full fre- 
quency dependence of the real and imaginary parts of the 
Matsubara-axis self-energy at hole doping x = 0.058 for 
the sectors (tt/2, tt/2), (0,7r), and (tt,tt/2). We can see 
immediately that for all sectors the criterion for Fermi- 
liquid behavior is violated. However, for the (tt/2, tt/2) 
and (0,7r) sectors, both the real and imaginary parts of 
the self-energy behave sufficiently smoothly that a rea- 
sonable extrapolation to uj = is possible with, in par- 
ticular, answers which are reasonably consistent whether 



the extrapolation is performed with the two lowest Mat- 
subara points or the three lowest. On the other hand, in 
the case of the (tt/2,tt) sector both the real and imagi- 
nary parts of the self-energy vary rapidly near lu = 0. 

The differences in Matsubara-axis self-energy are re- 
flected in differences in the maximum entropy analyti- 
cally continued self energies shown in Fig. 1151 (Other 
analytically continued results are shown In Fig. [TO]) We 
see that values of ImSjf(w = 0) for the K = (tt/2, tt/2) 
and (0,7r) sectors are approximately —1 and —2, very 
consistent with the values obtained from our poor-man's 
analytical continuation, whereas the ImS in the other 
sectors in fact displays a gap, which is not consistent 
with an interpretation of the extrapolated imaginary axis 
quantity as a scattering rate. 




FIG. 15: ImS(oj) for 16 site cluster at f3t = 7.5 and doping 
x — 0.058 obtained from maximum entropy analytical contin- 
uation of self-energy data shown above. 

The differences in continued behavior can be related to 
the structure of the Matsubara-axis data at low frequen- 
cies. In particular, we find that if extrapolations using 
two or three Matsubara-axis data points or three data 
yield answers which are consistent with about 5% and if 
the three lowest Matsubara points are all below the max- 
imum in I ImS I then in the quasi- Fermi-liquid regime (as 
defined by T < ttT) our poor man's analytical continu- 
ation procedure yields £"(£;, 0) values which agree with 
analytical continuation to within about 20—40% and Z^'s 
which are about 40% larger than the analytical continu- 
ation values, but with the same doping dependence. The 
scattering rate data we discuss in the text fulfill these 
conditions. Larger differences between extrapolated and 
analytically continued results occur outside of the Fcrmi- 
liquid regime and for data which violate the conditions 
outlined above. 



1 M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998). 

2 P. W. Anderson, Science 235, 1196 (1987). 

3 W. F. Brinkman and T. M. Rice, Phys. Rev. B 2, 4302 



(1970). 

4 A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

5 W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 



15 



(1989). 

6 Z. M. Yusof, B. O. Wells, T. Valla, A. V. Fedorov, P. D. 
Johnson, Q. Li, C. Kendziora, S. Jian, and D. G. Hinks, 
Phys. Rev. Lett. 88, 167006 (2002). 

7 M. Abdel-Jawad, J. G. Analytis, L. Balicas, A. Carrington, 
J. P. H. Charmant, M. M. J. French, and N. E. Hussey, 
Phys. Rev. Lett. 99, 107002 (2007). 

8 M. J. French, J. G. Analytis, A. Carrington, L. Balicas, 
and N. E. Hussey, New J. Phys. 11, 055507 (2009). 

9 T. Valla, A. V. Fedorov, P. D. Johnson, B. O. Wells, S. L. 
Hulbert, Q. Li, G. D. Gu, and N. Koshizuka, Science 285, 
2110 (1999). 

10 L. B. Ioffe and A. J. Millis, Phys. Rev. B 58, 11631 (1998). 

11 A. G. Loeser, Z.-X. Shen, D. S. Dessau, D. S. Marshall, 
C. H. Park, P. Fournier, and A. Kapitulnik, Science 273, 
325 (1996). 

12 H. Ding, T. Yokoya, J. C. Campuzano, T. Takahashi, 
M. Randeria, M. R. Norman, T. Mochiku, K. Kadowaki, 
and J. Giapintzakis, Nature 382, 51 (1996). 

13 A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. 
Phys. 75, 473 (2003). 

14 M. Grilli and G. Kotliar, Phys. Rev. Lett. 64, 1170 (1990). 

15 M. Grilli, B. G. Kotliar, and A. J. Millis, Phys. Rev. B 42, 
329 (1990). 

16 X.-G. Wen and P. A. Lee, Phys. Rev. Lett. 76, 503 (1996). 

17 P. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 
78, 17 (2006). 

18 A. Kampf and J. R. Schrieffer, Phys. Rev. B 41, 6399 
(1990). 

19 Y. M. Vilk, Phys. Rev. B 55, 3870 (1997). 

20 A. Abanov, A. V. Chubukov, and J. Schmalian, Europhys. 
Lett. 55, 369 (2001). 

21 B. Kyung, V. Hankevych, A.-M. Dare, and A.-M. S. Trem- 
blay, Phys. Rev. Lett. 93, 147004 (2004). 

22 C. Castellani, C. Di Castro, and M. Grilli, Phys. Rev. Lett. 
75, 4650 (1995). 

23 M. H. Hettler, A. N. Tahvildar-Zadeh, M. Jarrell, et al., 
Phys. Rev. B 58, R7475 (1998). 

24 M. H. Hettler, M. Mukherjee, M. Jarrell, and H. R. Krish- 
namurthy, Phys. Rev. B 61, 12739 (2000). 

25 A. I. Lichtenstein and M. I. Katsnelson, Phys. Rev. B 62, 
R9283 (2000). 

26 G. Kotliar, S. Y. Savrasov, G. Palsson, and G. Biroli, Phys. 
Rev. Lett. 87, 186401 (2001). 

27 T. Maier, M. Jarrell, T. Pruschke, and M. H. Hettler, Rev. 
Mod. Phys. 77, 1027 (2005). 

28 G. Kotliar, S. Y. Savrasov, K. Haule, et al., Rev. Mod. 
Phys. 78, 865 (2006). 

29 A. M. S. Tremblay, B. Kyung, and D. Senechal, Low Temp. 
Phys. 32, 424 (2006). 

30 C. Huscroft, M. Jarrell, T. Maier, S. Moukouri, and A. N. 
Tahvildarzadeh, Phys. Rev. Lett. 86, 139 (2001). 

31 A. Macridin, M. Jarrell, T. Maier, P. R. C. Kent, and 
E. D'Azevedo, Phys. Rev. Lett. 97, 036401 (2006). 

32 O. Parcollet, G. Biroli, and G. Kotliar, Phys. Rev. Lett. 
92, 226402 (2004). 

33 M. Civelli, M. Capone, S. S. Kancharla, O. Parcollet, and 
G. Kotliar, Phys. Rev. Lett. 95, 106402 (2005). 

34 C. Berthod, T. Giamarchi, S. Biermann, and A. Georges, 
Phys. Rev. Lett. 97, 136401 (2006). 

35 B. Kyung, S. S. Kancharla, D. Senechal, A.-M. S. Trem- 
blay, M. Civelli, and G. Kotliar, Phys. Rev. B 73, 165114 
(2006). 

36 T. D. Stanescu and G. Kotliar, Phys. Rev. B 74, 125110 



(2006). 

37 M. Ferrero, P. S. Cornaglia, L. D. Leo, O. Parcollet, 
G. Kotliar, and A. Georges, Europhys. Lett. 85, 57009 
(2009). 

38 M. Ferrero, P. S. Cornaglia, L. De Leo, O. Parcollet, 
G. Kotliar, and A. Georges, Phys. Rev. B 80, 064501 
(2009). 

39 M. Ferrero, O. Parcollet, A. Georges, G. Kotliar, and D. N. 
Basov, Phys. Rev. B 82, 054502 (2010). 

40 P. Werner, E. Gull, O. Parcollet, and A. J. Millis, Phys. 
Rev. B 80, 045120 (2009). 

41 E. Gull, O. Parcollet, P. Werner, and A. J. Millis, Phys. 
Rev. B 80, 245102 (2009). 

42 N. Lin, E. Gull, and A. J. Millis, Phys. Rev. B 82, 045104 
(2010). 

43 M. Jarrell, T. Maier, C. Huscroft, and S. Moukouri, Phys. 
Rev. B 64, 195130 (2001). 

44 S. Moukouri and M. Jarrell, Phys. Rev. Lett. 87, 167010 
(2001). 

45 T. A. Maier, M. Jarrell, T. C. Schulthess, P. R. C. Kent, 
and J. B. White, Phys. Rev. Lett. 95, 237001 (2005). 

46 P. R. C. Kent, M. Jarrell, T. A. Maier, and T. Pruschke, 
Phys. Rev. B 72, 060411 (2005). 

47 E. Kozik, K. Van Houcke, E. Gull, L. Pollet, N. Prokof'ev, 
B. Svistunov, and M. Troyer, Europhys. Lett. 90, 10004 
(2010). 

48 E. Gull, P. Werner, O. Parcollet, and M. Troyer, Europhys. 
Lett. 82, 57003 (2008). 

49 G. Alvarez, M. S. Summers, D. E. Maxwell, M. Eisenbach, 
J. S. Meredith, J. M. Larkin, J. Levesque, T. A. Maier, 
P. R. C. Kent, E. F. D'Azevedo, et al., in SC '08: Proceed- 
ings of the 2008 ACM/IEEE conference on Supercomput- 
ing (IEEE Press, Piscataway, NJ, USA, 2008), pp. 1-10, 
ISBN 978-1-4244-2835-9. 

50 O. K. Andersen, O. Jepsen, A. I. Liechtenstein, and I. I. 
Mazin, Phys. Rev. B 49, 4145 (1994). 

51 A. Comanac, L. de' Medici, M. Capone, and A. J. Millis, 
Nat Phys 4, 287 (2008). 

52 N. Trivedi and M. Randeria, Phys. Rev. Lett. 75, 312 
(1995). 

53 E. Gull, P. Werner, X. Wang, M. Troyer, and A. J. Millis, 
Europhys. Lett. 84, 37009 (2008). 

54 H. Park, K. Haule, and G. Kotliar, Phys. Rev. Lett. 101, 
186403 (2008). 

55 T. A. Maier, T. Pruschke, and M. Jarrell, Phys. Rev. B 
66, 075102 (2002). 

56 T. Maier, T. Pruschke, and M. Jarrell, Physica C: Super- 
conductivity 388-389, 94 (2003), ISSN 0921-4534, pro- 
ceedings of the 23rd International Conference on Low Tem- 
perature Physics (LT23). 

57 J. M. Luttinger, Phys. Rev. 119, 1153 (1960). 

58 M. Civelli, M. Capone, A. Georges, K. Haule, O. Parcollet, 
T. D. Stanescu, and G. Kotliar, Phys. Rev. Lett. 100, 
046402 (2008). 

59 K. Haule and G. Kotliar, Phys. Rev. B 76, 104509 (2007). 

60 I. M. Vishik, W. S. Lee, F. Schmitt, B. Moritz, 
T. Sasagawa, S. Uchida, K. Fujita, S. Ishida, C. Zhang, 
T. P. Devereaux, et al., Phys. Rev. Lett. 104, 207002 
(2010). 

61 A. Liebsch and N.-H. Tong, Phys. Rev. B 80, 165126 
(2009). 

62 G. Sordi, K. Haule, and A.-M. S. Tremblay, Phys. Rev. 
Lett. 104, 226402 (2010). 

63 S. Sakai, Y. Motome, and M. Imada, Phys. Rev. Lett. 102, 



16 

056404 (2009). 65 X. Wang, E. Gull, L. de' Medici, M. Capone, and A. J. 

A. Albuquerque, F. Alet, P. Corboz, et al., J. Magn. Magn. Millis, Phys. Rev. B 80, 045101 (2009). 

Mater. 310, 1187 (2007). 



